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1. Introduction 

The aim of this note is to present a modification of the standard Coordinate Bethe Ansatz so as to deal 
with open spin chains with non-diagonal boundary matrices. We applied this new method successfully 
to two models |T,T|, the XXX and XXZ open spin chains, but it clearly should work for other integrable 
models. To be as simple as possible, we will mainly stick to the XXX model. It will allow us to present in 
a pedagogical way the generalization we use. The XXZ case is then studied focusing on the differences 
with the previous case. 

The plan of the article is the following. We first introduce in section [2] the XXX model with periodic 
boundary conditions. It allows us to present the standard Coordinate Bethe Ansatz and to fix notations. 
Then, in section [3] we deal with the open XXX model, with one triangular boundary matrix. This 
case cannot be dealt with the standard Coordinate Bethe Ansatz, and we use a generalization of it. In 
section H] we study the XXZ model with non-diagonal boundary matrices: our method allows to get the 
eigenf unctions for boundary matrices obeying some constraints, some of them being known, the others 
new. In section HI we present another ansatz we recently introduced, the Matrix Coordinate Bethe Ansatz 
IJl- Together, the generalized Coordinate Bethe Ansatz and this new Ansatz provide the complete set of 
eigenvectors and eigenvalues for the XXZ open chain. We conclude in section [6] 



2. XXX model with periodic boundary conditions 

The XXX spin chain |4| is one of the most studied integrable models. To fix the notations, we start with 
the XXX model with periodic boundary conditions. It describes the interaction of spins ^ on a Id lattice 
(of L sites), with Hamiltonian 
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where I is the identity and Pe/+i the permutation operator acting on sites {£, £ + I): 
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P{u(g)v) =v(^u. (2) 



Here and below, we use auxiliary space notation: indices indicate on which sites of the chain operators 
act non trivially. For instance 

P34 = I(S)I0 P e End{-H) . (3) 

The interaction is a nearest neighbours interaction, and we set L + 1 = 1 (periodic boundary conditions). 
The Hamiltonian acts on an Hilbert space: Ti = (C^)^^, whose states take the form: 

1 1) ® 1 1) ® u> ® • • • 1 1) ® u> = I m • • • n) (4) 



2.1. Integrability 

It is well-known that the Hamiltonian is integrable and can be obtained from a so-called transfer matrix, 
see e.g. fS ] and references therein. Without going into details, we just remind the steps to get it. One has 
first to consider the Hamiltonian H = H + LI. It is clear that H and H have the same eigenfunctions. 
Then, one defines the transfer matrix 

t(A) = tro(i?oi(A)i?o2(A)---i?oL(A)) with Rm{\) = + Pu (5) 

where A is the spectral parameter, and P is defined in Q. The R matrix Rki{^) obeys the celebrated 
Yang-Baxter equation 

-Rl2(Al — A2) Rvi{\i — A3) i?23(A2 — A3) = R2'i{^2 — A3) i?l3(Ai — A3) i?l2(Ai — A2) • (6) 

From this relation, it is easy to see that [t(A) , t(A')] = 0, VA, A', so that upon expansion in A, t{\) 
generates L commuting independent charges. Since H = lni(A)|^_g, these charges are conserved, 
which proves the integrability of the model associated to H and H. 



2.2. Coordinate Bethe Ansatz 

We are looking for Hamiltonian eigenfunctions H <^ = E^. The solution to this problem using the 
Coordinate Bethe Ansatz has been known for a long time [6]. The starting point is a reference state that 
is a (zero energy) eigenstate: | t • • • t) = 0- 

Accordingly to this reference state, one can define a general state of the Hilbert space: 

|xi,...,x^) = |t...t4t---t4t t 4 t-.-t) € (C2)^^ (7) 

XI X2 Xm 

where xi < X2 < • • • < Xm are the positions of the m spins down of the state. Note that m is a quantum 
number: it corresponds to the operator ^ — where = X^^^^ s| is the z-component of the spin 
operator (we remind L is the number of sites). 

Then, the Coordinate Bethe Ansatz tlJ is a sort of plane-waves decomposition with respect to the 
above basis: 

E E 4'"^e*''^-"|2;i,---,^m). (8) 
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kj, j = 1, m are the plane wave momenta, Sm is the symmetric group {su{m) Weyl group), generated 
by transpositions aj, j = 1, . . . ,m — 1 that exchange kj and /cj+i, and 

^9 ~ (^g(l)' • • • T^g{m)) • (9) 

The coefficients A^g^^ are complex numbers to be determined such that 

H^m = E^<^m. (10) 

We project equation (ITOl ) on the different independent vectors . . . ,Xm) to get constraints on the 
coefficients aIj^\ Due to the form of H, it is enough to consider three cases only: 

• all the Xj's are far away one from each other (1 + Xj < Xj+i, Vj) and are not on the boundary sites 
1 and L. This case will be called generic. 

• Xj + 1 = Xj+i for one given j, 

• xi = 1 or Xm = L (periodicity condition L + 1 = 1). 

By linearity, more complicated cases just appear as superposition of 'simple' ones. 

These three projections lead to equations that one needs to solve. We do not reproduce the calculations 
here, but just give the solutions. 

> Calculation of the energy: projection on . . . , Xn) generic 

E^ = y" A(e*'=0 where A(u) =u + --2= ~ ^ . (11) 

i=i 

> Scattering matrix: projection on . . . , xj, Xj+i = 1 + Xj, . . . , Xm) 
It provides the scattering matrix between pseudo-excitations. 

S{u,v) = (13) 

2u — uv — 1 

Eq. (fT2l) allows to express all the coefficients ^g*"^ in term of a single one, say A^^\ where id is the 
identity in Sm- 

> Bethe equations: projection on \xi . . . , Xm~i, L) 

This last constraint consists in the quantization of the pseudo-excitation momenta kj since the system is 
in a finite volume. 

m 

-Q g^^ik, ^ ^ik, ^ ^ ^iLk, f^^j. (14) 
£=1 



With these three projections, one gets all the relevant physical information for the model, and obtains 
the eigenfunctions of the Hamiltonian. 



3. XXX model with boundaries 

The open XXX model has Hamiltonian 



L-1 L-1 

H = + Hbuik + Bl where Hbuik = X] = ^ " (15) 

with boundary matrices B^ that we choose of the form 

This Hamiltonian describes the interaction of spins (up or down) among themselves, and with two 
boundaries described by the matrices B^. These matrices preserve integrability of the model (see below). 

Let us stress that the boundary matrices corresponds to a new case, the B^ matrix being triangular 
(when fi 7^ 0). This means that the left boundary can now flip the spin 4, to This has drastic 
consequences as we shall see. 



3.1. Gauge transformations 

Obviously, any Hamiltonian H' related to i7 by a gauge transformation, 

i7' = [/ [/ (g, . . . [/ i7 [/-I (g, [/-I . . . [7-1 (17) 

^ V ' V ' 

L L 

will have the same spectrum as H, and their eigenfunctions will be related in an obvious way. 

Since is invariant under these gauge transformations, our approach is vahd for any boundary 

matrices deduced from B^ by a gauge transformation = U B^ U~^. 

In particular, for the case 



U = t\\ I ] (18) 



one gets lower triangular matrices 



A- = ( ^ f „, M and A- = ( , ' . M. (19) 

On the contrary, when a = /3 and 5 7^ 7, it is not possible to find a gauge transformation that 
diagonalizes B^ while keeping B diagonal, since the two matrices do not commute. 

It should be also clear that the same treatment can be done when B^ is diagonal and B is triangular. 



3.2. Integrability and connection with reflection algebra 

For those famihar with the so-called reflection equationlH |9l, let us note that the matrices B^ do not 
obey this equation. However, they are connected to such 'reflection matrices' in the following way. 
One first has to make them traceless, using the identity matrix (which does not change the form of the 
eigenfunctions) and consider the Hamiltonian 

H = H- ^{trB-) h - ^itrB+) U + (L - 1)1 = H + {L - 1 - ^{trB~) - ^{trB+)) I (20) 

where the indices 1 and L (that are in fact not relevant when considering the identity matrix) are explicited 
for obvious reasons. Then, the two new matrices B^ = B^ — ^{trB^) I are connected to 'reflection 
matrices' through 

K^{X)=I + XB^ (21) 



where A is the spectral parameter. These 'reflection matrices' obey 

i?i2(Ai - A2) K^iXi) i?i2(Ai + A2) K^{X2) = K^{\2) Ri2{\i + A2) i?i2(Ai - A2) Kf{\i) (22) 

where i2i2(A) has been defined in section |2?T] Again, from the Yang-Baxter equation ^ and the 
reflection equation (l22l ). one proves that the transfer matrix [91 

t{X) = tro(^K+{X) Roi{X) • • • Rol{X) K-{X) Rol{-X) ■ ■ ■ iioi(-A)) (23) 

obeys [t{X) , t{X')] = 0, VA, A'. Since one has H = ■^t{X)\^_Q, this proves the integrabiUty of the 
model associated to H. 

3.3. Generalized Coordinate Bethe Ansatz 

Again, as for the periodic case, the starting point is a reference state: | t • • • t) = (o^ + 7) I t • • • t)- 
Note that even when /U ^ this state is a reference state, while | | . . . |) is not anymore. 

> When /X = (diagonal boundaries) the boundaries do not modify the spin (no flip) and one can use 
the "usual" Coordinate Bethe Ansatz: 

E E 4'"^e^*'^"l^i,---,^m), (24) 

x\<---<Xm geBCm 

BCrn is the Bm Weyl group, generated by the symmetric group Sm and the reflection Ri exchanging ki 
and —ki. 

> When /U / 0, one has to modify the Ansatz 

n 

^n=Y: E E 4"-'"^ 1^-+!' ■■■'^n), (25) 

m=0 Xm+i<---<Xn geGm 

where G„ = BCJBC^ and fcf ^ • = Ej=,n+i 

Let us stress that in this model the number of pseudo-excitations i is not conserved, although the 
model is still integrable. 

The coefficients aIP'™'^ are all determined (but one) by 

H<fn = En^n. (26) 

We project this equation on states \x) with: 

• (xi,X2, ...,x„) generic (1 + xj < xj+i, Vj) 

• + 1 = Xj+i for some j 

• — L 

• Xi = \ 

• ■■■,Xn) generic (m > 0) 

As in the periodic case, these projections lead to equations that have to be solved. We do not reproduce 
the (rather lengthy) calculations, and only give the results. We refer to [2|| for details on the calculation. 
The results contain the relevant physical information of the model. 

> Calculation of the energy: projection on . . . , x„) generic 

^„ = a + 7 + V A(e^''0 where Xiu) = u + --2= ^'^~ ' . {11) 
^-^ u u 

i=i 

Note that it has a "bulk part" similar to the periodic case, and a boundary contribution that is independent 
of /i. 



> Scattering matrix: projection on |xi, . . . , Xj,Xj+i = 1 + Xj, . . . , Xn) 

4^-0) = sfe^^^PW^e^Wi)) 4^^'°) where Siu, v) = - ~ ~ ^ (28) 
^ ^ \ ) ^ 2u — uv — 1 

It is similar to the periodic case since the boundaries are not involved in this process. 

> Reflection coefficient for the left boundary: projection on [1, Xm+i • • • , Xn) 

A^-f = R(e^^.i) ^(«.o) where R{z) = -z^ ^'1 + ^-^ ^ !±(iM ^ ^^9) 
3^ ^ 1 — z + p — a Tj^{z) 

This equation is specific to open case (i.e. in presence of a boundary), but is valid whatever the boundary 
matrices are (diagonal or not). 

> Transmission coefficient: projection on \xm+i . . . , x„) 

J^{n,m) _ |-gjfc<,(i) ^ e^'^s*™' ) 4"'"^"^) (31) 

r(™)(zi,...,z^) = ^ (32) 

a(^i,^;2)=« • (33) 

ZiZ2 - 1 

This equation is specific to the case of triangular boundary matrices. It is new with respect to the case 
of diagonal boundary matrices. It relates the coefficients Ag"'™^ with different m's: it shows that the 
number of spins down cannot be conserved. In other words we have a system where the number of 
"pseudo-particles" (spin down) is not conserved, while the model is still integrable. 

> Bethe equations: projection on\x\ . . . ^ Xn-i, L) 

TT S{e'''' , e'^^ ) S{e-'^^ , e'^') = e^'^^^ _Pj M . l<3<n (34) 

<¥) 

r.{z) = ^{l-z + 6-j). (35) 

Z ~v ± 



Note that, as in eq. (I27I ). the Bethe equations do not depend on /i. This proves that the eigenvalues 
are the same as the ones of the model associated to diagonal boundary matrices. This correspondence 
ensures that the eigenvalues are real although the Hamiltonian is not Hermitian. 

Let us stress that although the eigenvalues do not depend on fi, the eigenvectors do. Hence the 
physical properties of the model are different. 

4. Generalization to XXZ model with boundaries 

The resolution of XXZ model with non-diagonal matrices shares the same ideas but with extra new 
features. Thus, we will not describe the approach in details and rather focus on these extra features, 
referring to [1] for details. However, to stick to the presentation done for XXX model, we will use XXZ 
notation, instead of the ASEP one used in U. The explicit form of the transformation relating the two 
notations can be found in e.g. llTOl . 



The Hamiltonian we consider has the form 



A = i(Q + Q-i) and (, = 1(Q - Q"!) (37) 
B ^ ( "..s -^r' ) \ -^T' ^ (38) 



where a are the usual Pauh matrices and // is a free parameter. 



4.1. Constraints between the (non-diagonal) boundary matrices 



To solve the XXX model, we considered one triangular and one diagonal boundary matrix (up to gauge 
transformations). In the same way, when one deals with the XXZ model, we need to consider special 
(non-diagonal) matrices of the form (1381 ). However, in addition to this special form, they have to obey 
some constraint relations: 

JJ (c,(q,7)c,/(/3,5) -Q^-^-"e-") =0 with c+{u,v) = - and c_(u,v) = 1 (39) 

or 

J] (2'e(a,7)2'e'(/3,<5)-g-"e^) = (40) 

~ , . -Q + v-u±^[Q-^-Q + v-uY + Auv 
with c±{u,v) = — ^ (41) 

where the integer n corresponds to the eigenfunction on which they act. Indeed, the first choice of 
constraints ( [39l ) has to be related to the original approach [11] that allowed to compute eigenvalues for 
XXZ model with non-diagonal boundary matrices using fusion relations. The second choice (|40l)-(l4T]) 
corresponds to new constraints. In both cases, we computed the eigenvalues and the eigenfunctions of 
the corresponding model. Below we present some results for the first choice of constraints, the complete 
treatment being done in HI. 



4.2. Basis vectors depend on which site they are 



The usual spin up and spin down vectors used in the XXX models have now to be replaced by the 
following vectors 

Remark that they depend on the site I where they stand and also of extra parameters Uk,dk- This 
site dependence has to be related to the local gauge transformations lfT2l that are used to construct the 
Algebraic Bethe Ansatz for XXZ model with non-diagonal matrices. 
Then, a generic Hilbert space vector is defined by 

l^m+l ) • • • ) ^n) 

and the parameters Um+i, Un and dm+i, dn are fixed by the generahzed Coordinate Bethe Ansatz. 
In particular they obey the relations ug^i = Q~^U£ and = Q~^di. 



4.3. Telescoping terms appear 

When local Hamiltonians /i^^+i act on generic vectors ( |43] ). they make appear new vectors \t) 
that are not of the form i 



hi2\d)®\u) = \d) ^\u) -Q\u) ®\d) + {Q -Q-'^)\d) ®\t) (44) 
hi2\u)(^\d) = \u)(^\d)-Q-^\d)(^\u)-{Q-Q-'^)\d)®\t) (45) 

hi2\d)(^\d) = {Q-Q-^)[\d)(^\t) -\t)®\d)'^ (46) 

hi2\u)®\u) = (47) 

In view of these relations, one could be tempted to think that the basis ( |43] ) is not suited for the study of 
the XXZ Hamiltonian. However, this is only true for local Hamiltonians On the general Hamil- 

tonian (l36l ). these new vectors appear with alternating signs, so that their only contribution to the total 
Hamiltonian is on the first and last site, where they are used to diagonalize the boundary matrices. 

Apart from these three modifications, the generalized Coordinate Bethe Ansatz for the XXZ model fol- 
lows the same steps as for the XXX one. 



5. Completeness and Matrix Coordinate Bethe Ansatz 

In the case of XXX model, it has been shown that the Coordinate Bethe Ansatz provides the complete set 
of eigenvectors for periodic boundary conditions |fT3l . It is also believed to be complete for open diagonal 
boundary conditions. For triangular boundary matrix, the spectrum being the same as for diagonal ones, 
the set should be complete too |2|. 

For XXZ model, it is known that the generalized Coordinate Bethe Ansatz do not provide all the 
eigenvectors. For instance, by numerical investigations [14], it has been established that the whole 
spectrum is given by two different types of Bethe equations. The present method provides the eigenvalues 
and eigenvectors corresponding to only one type of Bethe equations. Remark that since the Hamiltonian 
is not Hermitian, the right and left eigenvectors are different. Thus, using the Coordinate Bethe Ansatz 
on left vectors leads to another set of eigenvalues. Together, the "right" and "left" eigenvalues generate 
the complete spectrum. However, in this way, one constructs only "half of the eigenvectors for each 
"side" (right or left vectors). Although this generalized Coordinate Bethe Ansatz is not enough to obtain 
all the vectors, it has the advantage of giving an interpretation of the number n entering in the constraint: 
it is the maximal number of pseudo-excitations in the Ansatz. 

To get a complete set of eigenvectors, one needs to use another ansatz: in few cases, one can use the so- 
called Matrix ansatz ifTSi (used in Statistical Physics), but in general it is not sufficient since it provides 
only one eigenvector. 

In 1 3 1, we developed a Matrix Coordinate Bethe Ansatz, that is a mixing of generalized coordinate 
Bethe ansatz and of Matrix ansatz. More precisely, it is a non-commutative generalized coordinate Bethe 
ansatz, where the entries in (1431 ) now belong to an algebra (very closed to the one introduced in [15 ]) 
acting in an additional auxiliary space. In this framework, the Matrix ansatz eigenvector appears as a 
new vacuum on which we build this non-commutative generalized coordinate Bethe ansatz. Numerical 
studies indicate the spectrum is then complete. For more details, we refer to the recent work 



6. Conclusion 

We have shown a generalization of the Coordinate Bethe Ansatz that allows to take into account the case 
of non-diagonal boundary matrices. The Ansatz has been applied to XXX and XXZ open spin chains, 
but it should also work on different integrable models. In the case of XXZ model, the Ansatz allows to 
recover and generalize the constraints found with different methods. However, the case of fully general 
boundary matrices remains to be done. 

In the case presented here, completeness of the Ansatz is ensured by the introduction of another 
Ansatz, the Matrix Coordinate Bethe Ansatz. A synthetical presentation of both Ansatzen is also lacking 
for the moment. 

Finally, let us stress that these Ansatzen can also be applied to open spin chains built on algebras of 
higher rank. 

Works are in progress on these subjects. 

Bibliography 

[1] N. Crampe, E. Ragoucy and D. Simon, JSTAT 1011 (2010) P11038 and arXiv: 1009.4119. 

[2] N. Crampe and E. Ragoucy, arXiv: 1 105.41 19. 

[3] N. Crampe, E. Ragoucy and D. Simon, J. Phys. A44 (2011) 405003 and larXiv: 1 106.32641 

[4] W. Heisenberg, Zeitschrift fur Physik 49 (1928) 619. 

[5] L. Faddeev, in Symetries Quantiques, Les Houches summerschool proceedings 64, Eds A. Connes, K. Gawedzki and J. 
Zinn-Justin, North-Holland 1998 and hep-th/ 9 6051 87. 

[6] M. Gaudin, Phys. Rev. A4 (1971) 386. 

[7] H. Bethe, Zeitschrift fur Physik 71 (1931) 205. 

[8] I.V. Cherednik, Theor. Math. Phys. 61 (1984) 977. 

[9] E.K. Sklyanin, J. Phys. A21 (1988) 2375. 

[10] J. de Gier and R Essler, J. Stat. Mech. (2006) P1201 1 and arXiv:cond -mat/0609645 1 

[11] R. Nepomechie, J.Phys. A34 (2001) 9993 and arXiv:hep-th/01 10081. 

[12] J. Cao, H. Lin, K. Shi and Y. Wang, Nucl. Phys. B663 (2003) 487 and arXiv:cond-mat/0212163. 

[13] A. N. Kirillov, Zap. Nauch. Sem. LOMI 131 (1984) 88. 

[14] R.I. Nepomechie and R Ravanini, J. Phys. A36 (2003) 1 1391 and arXiv:hep- th/0307095| 

[15] B. Derrida, M. Evans, V. Hakim and V. Pasquier, J. Phys. A26 (1993) 1493. 



